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ABSTRACT 

A  study  of  the  modifications  produced  by  high  power,  high  frequency  (HF) 
waves  in  the  F-region  of  a  quiet  (low  solar  cycle)  polar  ionosphere  is  pre¬ 
sented.  In  the  polar  geometry,  maximum  Ohmic  heating  occurs  producing  signi¬ 
ficant  changes  to  the  zero  order  electron -temperature  and  density  profiles. 

The  temporal  and  spatial  aspects  of  these  changes  are  calculated  and  their 
effects  on  relevant  ionospheric  parameters,  such  as  the  d.c.  conductivity  ten¬ 
sor,  are  found.  An  Interesting  finding  is  the  generation  of  density  cavities 
at  the  reflection  point  of  the  HF  waves.  A  survey  of  the  dependence  of  the 
various  effects  on  HF  pump  parameters  is  made. 


I.  INTRODUCTION 


Presently,  there  is  renewed  interest  in  the  problem  of  ionospheric  mod¬ 
ification  by  high  power  HF  waves.  To  a  large  extent  this  activity  has  been 
stimulated  by  the  recent  completion  of  versatile  experimental  facilities  near 
the  polar  region,  namely,  HEATER  in  Norway  [Stubbe  et  al.,  1982]  and  HIPAS  in 
Alaska  [Hong  et  al. ,  1981].  Although  the  ionospheric  modification  problem  has 
been  extensively  investigated  over  the  past  decade  (see  for  reference,  Vol.  9_, 
issue  No.  11,  Radio  Science  [1974];  Vol.  44,  issue  No.  12,  Journal  of  Atmos¬ 
pheric  and  Terrestrial  Physics  [1982]),  the  majority  of  theoretical  and  nu¬ 
merical  investigations  have  been  primarily  directed  to  understanding  processes 
of  relevance  to  the  mid-latitude  ionosphere  [Melts  and  Lelevier,  1970;  Melts 
et  al. ,  1974] .  This  is  essentially  because  of  the  location  of  the  earlier  ex¬ 
perimental  facilities  in  Platteville,  Colorado  and  Arecibo,  Puerto  Rico.  Al¬ 
though  some  of  the  physical  processes  at  work  in  mid-latitude  ionospheric  mod¬ 
ification  are  also  encountered  in  the  polar  region,  the  latter  environment  has 
certain  unique  features  which  demand  specialized  attention.  With  this  per¬ 
spective  in  mind,  the  present  study  alms  to  establish  a  background  for  the 
simplest  study  of  ionospheric  modification  in  the  polar  F-reglon,  from  which 
more  sophisticated  models  can  develop.  It  is  hoped  that  this  study  may  serve 
as  a  reference  against  which  some  of  the  ongoing  experiments  can  be  compared 
as  well  as  a  means  of  assessing  the  potential  of  future  field  studies. 

Among  the  various  unique  features  encountered  in  the  polar  ionosphere, 
the  fact  that  the  magnetic  field  of  the  earth  makes  a  small  angle  with  the 
plasma  density  gradient  has  a  profound  imact  on  the  heating  produced  by  HF 
waves.  The  reason  is  that  the  heating  source  (i.e. ,  the  HF  beam)  can  be 
nearly  colinear  with  the  direction  of  the  heat  flow  (i.e.,  parallel  to  the 


■agnetic  field  lines) ,  resulting  in  enhanced  heating  along  a  magnetic  flux 
tube.  This  situation  is  to  be  contrasted  with  a  mid-latitude  or  equatorial 
heating  geometry  in  which  only  a  small  portion  of  a  flux  tube  is  exposed  to 
the  HF  beam.  The  results  of  the  present  Investigation  aim  to  emphasize  the 
role  of  the  collnear  heating  by  concentrating  on  the  case  7Neo  x  Bn  -  0,  where 
Neo  is  the  zero  order  plasma  density,  and  Bq  the  earth's  magnetic  field.  Some 
interesting  features  of  the  colinear  heating  geometry  have  already  been  exam¬ 
ined  in  the  context  of  simulating  microwave  beam  heating  from  power  station 
satellites  [Perkins  and  Roble,  1978].  A  technical  simplif ication  which  re¬ 
sults  from  the  consideration  of  the  idealized  VNeo  *  B0  ■  0  geometry  is  that 
the  spatial  dependence  of  the  HF  electric  field  can  be  well  represented  analy¬ 
tically,  thus  permitting  the  consideration  of  overdense  heating  scenarios, 
l.e.,  heating  below  the  critical  frequency  of  the  F-region  (u^).  Enhanced 
absorption  near  the  reflection  layer,  of  great  experimental  interest  [e.g. , 
Gordon  and  Carlson,  1974;  Mantas  et  al. ,  1981],  is  thus  accurately  calculated 
without  recourse  to  the  approximations  necessary  in  mid-latitude  studies 
[Thomson,  1970;  Meltz  et  al. ,  1974] .  Also,  since  in  this  geometry  linear  mode 
conversion  and  parametric  decay  instabilities  are  absent,  the  only  heating 
channel  is  purely  Ohmic  (colllsional) .  Therefore,  by  emphasizing  this  geome¬ 
try  one  obtains  an  upper  bound  on  the  magnitude  of  the  effects  resulting  from 
Ohmic  heating. 

The  present  spatial  and  temporal  study  of  transport  in  a  model  polar  ion¬ 
osphere  addresses  only  the  changes  produced  in  the  F-region  plasma.  The  com¬ 
plicating  effects  arising  from  the  chemistry  of  the  D  and  E-regions  are  best 
studied  separately  since  they  occur  on  different  time  scales,  and  are  not  sen¬ 
sitive  to  the  unique  VNeo  *  Bq  *  0  polar  geometry  because  of  the  low  degree  of 


magnetization  at  lower  heights.  In  addition,  many  of  the  experimental  effects 
presently  being  studied  take  place  in  the  F-region  [e.g. ,  Mantas  et  al. ,  1981; 
Stubbe  et  al. ,  1982] .  Hence  it  is  of  importance  to  assess  the  role  of  colli- 
sional  heating  in  this  region. 

In  describing  the  interaction  of  HF  heated  electrons  with  the  various 
neutral  components,  the  present  investigation  incorporates  the  substantial 
body  of  knowledge  available  in  this  area;  see  for  example  Schunk  and  Nagy 
[1980]  and  the  references  therein.  In  particular,  this  work  relies  heavily  on 
the  works  of  Banks  and  Kockarts  [1973]  and  Rees  and  Roble  [1975]  in  evaluating 
the  collision  cross  sections,  the  elastic  and  inelastic  cooling  losses,  and 
the  thermal  conductivity. 

The  present  investigation  uses  an  analytic  expression  for  the  electric 
field  of  the  HF  pump  valid  for  the  polar  F-region.  The  Ohmic  absorption  of 

i 

this  field  is  calculated,  as  well  as  the  resulting  absorption  coefficient. 

The  self-consistent  evolution  of  the  temperature  and  density  variations  is  nu¬ 
merically  calculated.  Their  feedback  on  relevant  ionospheric  parameters  is 
found.  A  survey  is  made  of  the  dependence  on  background  density  and  on  HF 
heater  parameters  such  as  power,  frequency,  polarization,  and  pulse  modulation 
(or  pulse  length). 

The  manuscript  is  organized  as  follows.  In  section  II,  the  four  parts  of 
the  model  are  described:  the  polar  F-region,  the  heat  transport  phenomena, 
the  pump  model,  and  the  method  of  solution.  Section  III  presents  the  results 
in  four  main  parts:  the  response  to  HF  heating,  power  dependence,  frequency 
dependence,  and  effect  of  HF  modulation  pulses.  A  discussion  of  the  results 
and  conclusions  are  presented  in  section  IV. 
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II.  HP  HEATING  MODEL 

A.  Polar  F-Region 

The  model  F-reglon  plasma  is  assumed  to  be  dominated  by  singly-charged 
atomic  oxygen  ions,  other  ionic  species  are  neglected.  The  major  role  of  the 
oxygen  ions  is  to  determine  the  ambipolar  rate  of  expansion  of  plasma  density 
resulting  from  HF  heating.  Of  course,  these  ions  also  determine  the  electron- 
ion  collision  rate  vej  which  is  the  dominant  heating  channel  in  the  upper  F- 
region.  Prior  to  irradiation  with  HF  waves,  the  ionosphere  is  assumed  to  be 
in  a  state  of  dynamic  equilibrium.  The  present  study  is  not  concerned  with 
determining  this  equilibrium  from  basic  solar  and  atmospheric  parameters,  as 
is  done  in  purely  ionospheric  modelling  [Schunk  and  Nagy,  1980],  but  rather  it 
is  assumed  that  the  zero  order  plasma  density  and  electron  temperature  pro¬ 
files  are  given.  The  profiles  selected  for  numerical  investigation  are  con¬ 
sistent  with  those  observed  experimentally  under  quiet  polar  conditions  [Bates 
and  Hunsucker,  1974].  The  density  and  temperature  profiles  selected  for  dis¬ 
cussion  in  this  paper  are  shown  in  Figure  la.  The  ion  temperature  is  taken 
to  be  equal  to  the  neutral  temperature  Tn.  At  the  lower  portion  of  the  F- 
region  the  electron  temperature  drops  to  the  neutral  temperature.  Throughout 
the  time  dependent  calculation,  both  the  ion  and  neutral  temperatures  are  kept 
constant  and  equal  to  each  other.  This  approximation  is  quite  reasonable 
since  the  HF  predominantly  heats  electrons,  due  to  their  much  lighter  mass. 

It  is  assumed  that  the  neutrals  in  the  F-region  are  adequately  described 
by  the  model  of  Banks  and  Kockarts  [1973].  The  neutral  components  are:  mol¬ 
ecular  nitrogen  (N2),  molecular  oxygen  (O2),  atomic  oxygen  (0),  atomic  hydro¬ 
gen  (H),  and  helium  (He).  The  altitude  concentrations  of  these  neutrals  are 
determined  from  the  Boltzmann  factor  with  a  constant  temperature  Tn  -  750°K. 
The  resulting  neutral  density  profiles  are  shown  in  Figure  lb.  From  the  plas- 


ma  density  and  temperature  profiles  of  Figure  la  and  the  neutral  profiles  of 
Figure  lb,  one  calculates  [Banks  and  Kockarts,  1973]  the  various  momentum 
transfer  collision  frequencies  entering  into  the  transport  description  of  HF 
heating,  namely,  electron-ion  collisions  ve^,  electron-neutral  collisions  ven, 
and  ion-neutral  collisions  v^n,  whose  profiles  are  shown  in  Figure  lc. 

B.  Heat  Transport 

Shortly  after  the  HF  pump  is  turned  on,  electrons  in  the  F-region  begin 
to  heat  up  due  to  momentum  transfer  collisions  with  neutrals  and  oxygen  ions. 
The  resulting  increase  in  electron  temperature  is  limited  by  two  different 
mechanisms  that  operate  simultaneously  but  with  varying  degree  of  importance 
depending  on  altitude.  The  two  mechanisms  are:  (1)  collisions  (due  to  both 
elastic  and  inelastic  processes)  with  neutrals  and  oxygen  ions,  and  (2)  field- 
aligned  electron  heat  conduction.  The  role  of  heat  conduction  across  the  mag¬ 
netic  field  is  relatively  unimportant  due  to  the  strong  electron  magnetization 
in  the  F-region  [Gurevich,  1978] ,  hence  it  is  neglected.  Furthermore,  heat 
convection  due  to  bulk  plasma  drifts  across  the  magnetic  field  becomes  impor¬ 
tant  for  time  scales  of  the  order  of  few  hundred  seconds  [Holway  et  al. , 

1978] .  Since  the  HF  pulse  durations  considered  in  this  study  are  shorter  than 
drift  times,  cross-field  heat  convection  is  ignored.  Although  the  neglect  of 
this  effect  is  justifiable  in  the  context  of  the  present  study,  its  role  be¬ 
comes  quite  important  when  attempting  to  describe  CW  conditions.  The  reason 
is  that  while  it  usually  takes  about  10  seconds  for  the  temperature  to  reach 
steady  state,  it  takes  several  minutes  for  the  density  to  diffuse.  Conse¬ 
quently,  the  plasma  drifts  away  from  the  HF  heated  volume  before  it  has  time 
to  reach  a  new  pressure  equilibrium.  The  result  of  this  effect  invalidates 


Che  density  calculations  of  prior  studies  based  on  steady  state  heating  of  the 
F-reglon  [Farley,  1963;  Gurevich,  1967;  Ignat'ev  et  al. ,  1977]. 


Since  electron  heat  conduction  takes  place  along  the  magnetic  field  and 
the  HF  pulse  considered  propagates  along  Bg  (for  the  case  VNeo  x  ^  ■  0)  the 
transport  equations  are  dominated  by  field  aligned  gradients.  The  weak  cross¬ 
field  gradients  arise  only  through  the  transverse  extent  of  the  HF  beam  (e.g. , 
the  Gaussian  cross  section  of  the  beam  spot).  Such  weak  variations  can  be 
treated  by  considering  each  field  line  individually  and  using  the  appropriate 
HF  electric  field  strength  on  each  line.  Although  this  procedure  is  relative¬ 
ly  straightforward  to  implement,  in  the  present  paper  we  report  only  on  the 
plasma  changes  taking  place  along  the  field  line  passing  through  the  center  of 
a  Gaussian  cross  section  HF  beam.  However,  the  effect  of  beam  divergence  is 
appropriately  taken  into  account  in  determining  the  HF  electric  fields,  as  de¬ 
rived  in  part  C  of  this  section. 

Using  a  coordinate  system  in  which  the  z  axis  points  along  the  direction 
of  increasing  height,  and  the  magnetic  field  of  the  earth  points  in  the  nega¬ 
tive  z  direction  (to  simulate  the  northern  hemisphere  situation),  the  one- 
dimensional  transport  equations  that  determine  the  self-consistent  evolution 
of  the  electron  density  Ne  and  temperature  Te  take  the  form 


25i  +  (2_  z)  •  (N  v  )  -  s  , 

at  3z  e  e  ° 


(i) 


N  v  -  _ ~~ _ 2_  [N  (T  +T  . )  ]  z  , 

e_e  <«even+mivin>  3z  e  e  1 


1  k  (N  22k  +  N  (v  #z)  lie}  +  IJ  kT  !_  (v  *z) 

2  ‘  e  3t  e  -e  3z  1  e  e  3z  -e  - 

2_[K  3Te_1  +  Q  +  Q  -  L 
3?«3T  o  e 


\  3T. 


(2) 


(3) 
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where  k  is  Boltzmann's  constant,  Vg  is  the  slow  electron  diffusion  velocity, 
mg  and  m^  are  the  electron  and  oxygen  ion  masses,  Ke  is  the  parallel  coeffi¬ 
cient  of  thermal  conductivity,  Q  is  the  HF  heat  source,  and  Le  represents  the 
rate  of  electron  cooling  due  to  both  elastic  and  inelastic  collisions.  The 
terms  SQ  and  Q0  account  for  the  steady  state  sources  that  bring  about  the 
background  equilibrium  profiles  in  the  absence  of  the  HF  pump.  The  quantities 
SQ  and  Q0  are  not  calculated  explicitly,  but  rather  they  can  be  obtained  by 
setting  3/3t  »  0  in  equations  (1)  -  (3)  in  which  the  profiles  of  Figure  la  are 
inserted.  We  assume  the  origin  z  ■  0  at  the  bottom  of  the  F-region,  which  is 
taken  to  represent  an  altitude  of  150  km. 

Equation  (2)  can  be  obtained  from  the  electron  and  ion  momentum  equations 
in  the  appropriate  slow  time  approximation  (i.e.,  neglecting  inertia)  and  for 
negligible  neutral  drift  velocity.  Equation  (3)  describes  the  local  conserva¬ 
tion  of  electron  energy  and  Includes  all  possible  channels  for  energy  gain  and 
loss.  The  second  and  third  terms  on  the  left  hand  side  of  (3)  represent  the 
effect  of  convection  and  pressure  flux,  -nd  are  found  to  be  negligible  in  the 
F-region.  However,  for  completeness  they  are  retained  in  the  numerical  calcu¬ 
lations. 

The  quantity  describing  the  rate  of  electron  energy  loss 
Le  =  Le(z,Ne,Te,T^,Tn)  contains  numerous  terms  not  worth  reproducing  here 
since  they  are  calculated  following  the  well  documented  works  of  Banks  and 
Kockarts  [1973],  Rees  and  Roble  [1975],  and  Schunk  and  Nagy  [1980].  However, 
it  is  worth  emphasizing  that  Le  includes  excitation  of  rotational  and  vibra¬ 
tional  levels  of  the  molecular  species  N2  and  O2,  excitation  of  atomic  oxygen 
fine  structure  0(^p)  at  the  ground  state,  and  excitation  of  atomic  oxygen  in  its 
first  excited  state  0(*D).  The  individual  terms  contributing  to  Le  are  shown 


in  Figure  2  for  F-region  parameters  prior  to  HF  heating.  The  importance  of 
the  inelastic  processes  in  the  lower  F-region  is  apparent,  specially  that  of 
the  excitation  of  the  0(3P)  level.  Electron-ion  collisions  become  the  major 
loss  mechanism  in  the  upper  F-region. 

The  coefficient  of  thermal  conductivity  used  in  (3)  for  a  partially- 
ionised  plasma  takes  the  form  [Banks  and  Kockarts,  1973,  equation  (22.116)] 

K  -  7-75»105Te5/2  .v  c.-l  „.c-l  -K-l  ,  W> 

6  1+3.22x104(T  2/N  )l  N  <Q  > 

e  e  n  n  n 

where  the  summation  in  the  denominator  is  over  all  the  neutral  gas  species  of 
concentration  Nn  and  average  electron-neutral  momentum  transfer  cross  section 
<Qn>.  The  Te5/2  dependence  is  associated  with  a  fully-ionised  plasma  and  it 
is  modified  by  the  denominator  term  to  account  for  the  effect  of  the  neutral 
gases. 

C.  HF  Pump 

The  source  term  Q  appearing  in  equation  (3)  is  related  to  the  HF  electric 
field  amplitude  Z  and  the  high  frequency  conductivity  tensor  a  through 

Q  “  <&'&>  *  y  ®  *  §]  ,  (5) 

where  the  brackets  <  >  imply  time  averaging,  £pthe  total  conduction  current, 
and  £  the  total  electric  field  whose  full  space-time  dependence  is  given  by 


£(z,t)  -  ^[E(z)e"iu,t+c.c.  ]  , 


(6) 
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in  which  u)  is  the  angular  frequency  of  the  HF  waves.  In  general,  the  total 
electric  field  can  be  described  in  terms  of  a  right  hand  circularly  polarized 
component  E+,  a  left  hand  circularly  polarized  component  E_,  and  a  parallel 
component  Ez.  In  the  present  geometry  (VNeo  x  =  0),  and  for  HF  beam  propa¬ 
gation  along  Bq,  only  the  E+  and  E_  components  are  relevant.  Accordingly,  the 
spatial  dependence  of  the  electric  field  amplitude  can  be  expressed  in  the 
form 


lS(z)  -  E+(z)  exp[i  / dz  k±(z)]  e+  ,  (7) 

where 

e±  -  —  (x  ±  iy)  ,  k±(z)  =  k Q)f  ±(z)  »  (8) 

/2 

where  kQ  *  ui/c,  and  k+  is  a  complex  quantity  in  general.  The  local  wavenumber 
k+(z)  can  itself  be  positive  or  negative  for  either  polarization  depending  on 
whether  the  wave  is  travelling  upwards  into  the  ionosphere  or  downwards  ifter 
reflection.  For  the  case  of  principal  interest  in  this  study,  i.e.,  overdense 
heating,  both  signs  of  <±  must  be  taken  into  account.  The  complex  index  of 
ref raction^/±(z)  -  n+(z)  +  ia+(z),  with  n+(z)  and  a+(z)  purely  real,  is  rela¬ 
ted  to  both  the  real  dielectric  and  HF  electric  conductivity  by 

z)  *  (n±(z)  +  ia±(z)]2  =  e±(z)  +  i  ^  o^z)  ,  (9) 

where  e±(z)  and  o+( z)  are  components  of  the  diagonal  dielectric  and  conductiv¬ 


ity  tensors,  respectively.  They  are  given  by 


(10a 
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the  parallel  components  of  these  tensors  are,  respectively. 


e  (z) 
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l  -  “p2(g> 


(u)2+V  2) 
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Cz)  ■  v«  "P  <*> 

4it  (uj2+v  2) 
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dob: 


where  <»»p(z)  is  the  local  electron  plasma  frequency,  toc  is  the  electron  gyro- 
frequency,  and  ve  is  the  total  electron  collision  frequency  for  momentum 
transfer,  i.e.,  ve  -  lven(Nn,  Te)  +  vgl  (Ng,  Te)].  For  F-region  propagation 


vg  <<  oi,  oic,  therefore,  away  from  the  cut-off  layer  a+2(z)  «  n±2(z), 


Accordingly , 


1/2 


n+(z)  51  [e+(z)  ]  ,  and  a+(z)  <* 


2n  o+(z) 


o)  n+(z) 
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The  WKB  description  of  (7),  although  quite  accurate  over  most  of  the 
travel  path  in  the  F-region,  needs  to  be  corrected  for  two  significant  ef¬ 
fects:  1)  beam  divergence  due  to  3-D  propagation;  2)  Alry-like  behavior  near 
the  turning  point  defined  by  e+(z=,z0)  =  0.  The  beam  divergence  correction  is 
introduced  by  assuming  a  conical  beam  spread,  as  illustrated  in  Figure  3,  and 
considering  the  reflection  layer  to  be  a  perfect  mirror  of  zero  width.  Con¬ 
servation  of  local  energy  flow  requires  that 


V  •  P  -  -  Q(z)  , 


(12 


where  is  the  Poynting  vector  whose  magnitude  is  P.  Integrating  (12)  over  a 
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volume  element  having  an  area  equal  to  the  HF  beam  cross  section  A(z)  yields 


P(z2)  A(z2)  -  P(zr)  A(z^ 


r 

-J dz  A(z)  Q(z) 
21 


(13) 


for  two  arbitrary  heights  z and  In  reality  the  beam  cross  section  A(z) 

corresponds  to  that  of  a  spherical  cap,  however  in  the  present  context  it  is 
approximated  by  the  area  of  a  disk,  as  is  appropriate  for  narrow  beams  pro¬ 
duced  by  high  gain  HF  antennas.  The  ratio  of  the  areas  is  equal  to 


A(z1)/A(z2)  -  (D  +  z1)2/(D  +  z 2> 2  , 


(14) 


where  D  is  the  distance  from  the  ground-based  antenna  to  the  bottom  of  the  F- 
region  at  150  km  altitude  (z  -  0),  as  indicated  in  Figure  3.  Incorporating 
the  beam  divergence  effect  into  equation  (7)  yields  an  expression  for  the 
electric  field  which  is  correct  in  the  region  0  <  z  <  zn,  where  zm  is  a  height 
slightly  below  the  turning  point  located  at  zG  as  illustrated  in  Figure  3, 


E±(z)  -  l'“) [e±(o)/e+(z)] 


1/4 


{ 


1/2 

[A(o) /A(z) ]  exp[ik 


,/dz.Jft 


♦(*')] 


1/2 

+  R[A(o)/A(2z0-z)] 


(15) 


Equation  (15)  contains  both  an  upward  travelling  wave  and  a  reflected  wave 
with  reflection  coefficient  R.  E(o)  is  the  electric  field  normalized  to  the 
HF  electric  field  delivered  from  a  ground-based  transmitter  to  the  bottom  of 
the  F-region,  i.e.. 


|E(o)  |  -  ( _ 2  )(^6°9WS)  V/m,  (16) 

l+nt(o)  d 
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where  U  is  the  average  power  radiated  in  watts,  G  is  the  antenna  gain,  D  is 
the  distance  from  the  ground  in  meters,  and  a  simple  dipole  radiator  is  as¬ 
sumed.  To  account  for  the  possibility  of  HF  beam  depletion  due  to  absorption 
in  the  D  and  E-regions  a  factor  0  <  0  <  1  is  introduced. 

The  correction  due  to  Airy-like  behavior  near  the  turning  point  is  imple¬ 
mented  by  having  a  separate  expression  for  the  electric  field  in  the  region 
zm  <  z  given  by 

E+(z)  -  2/x  (k0£)  ^  (e+(o)]  ^  E(o)  [A(o)/A(z)]  ^  Ai(c-ir) 

x  exp[ikQ  |  dzj(±(z)-iir/4] ,  (17) 

o 

where  C  ■  (k T  *  (kQ£)2/3  Im(jT±^(z) }»  and  it  is  assumed  that 

-1 

near  the  turning  point  e±(z)  -  -(z-z Q)/l,  where  l  =  [d(ln  Neo)/dz]z  is  the 

o 

local  density  scale  length  at  the  turning  point.  The  linear  density  approxi- 

2/3 

mation  is  not  a  limiting  factor  in  the  analysis  as  long  as  (kQi)  »  1, 
i.e.,  the  local  ionospheric  density  scale  length  (~  50  km)  is  much  larger  than 
the  width  of  the  main  lobe  of  the  Airy  pattern  (~  500  m).  The  linear  approxi¬ 
mation  would  only  be  invalidated  if  the  pump  frequency  and  the  critical  fre¬ 
quency  of  the  F-region  are  so  close,  that  the  cut-off  layer  and  the  critical 
frequency  layer  are  separated  by  a  distance  Az  equivalent  to  the  Airy  pattern's 
main  width,  i.e.,  Az  ~  i/(k0i)2/3  which  is  then  a  case  of  edge-diffraction  not 
of  interest  in  the  present  work.  However,  we  should  note  that  in  such  a  case 
a  parabolic  density  profile  at  the  cut-off  layer  is  a  more  appropriate  choice 
[Thomson,  1970;  Banos  and  Kelly,  1974], 

The  asymptotic  form  of  equation  (17)  matches  directly  to  the  WKB  solution 
of  equation  (15)  provided  the  reflection  coefficient  R  has  the  form  [Budden, 
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1961] 


R 


-i 


exp (-i2k 


(18) 


o 

Thus,  we  have  an  analytical  description  of  the  HF  electric  field  over  the  en¬ 
tire  F-region.  The  heat  source  Q  entering  in  equation  (3)  can  be  now  simply 
calculated  from  equations  (5),  (10),  and  (15)  or  (17). 

We  should  note  at  this  point  that  In  actuality  (kQi)  is  usually  so  large 
that  k0|z0  -  zm|  <<  kQi,  and  WKB  solutions  are  valid  practically  everywhere 
[Thomson,  1970].  The  full  wave  solution  is  therefore  confined  to  a  small  spa¬ 
tial  region  [Figure  3] .  However,  within  this  small  region,  large  amplitude 
gradients  appear  and  are  important  in  determining  the  shape  and  magnitude  of 
the  density  cavities  that  are  created  at  the  reflection  layer. 

The  total  amount  of  HF  power  absorption  pa  due  to  Ohmic  heating  for  the 
overdense  case  of  interest  here  is  given  by  equation  (13)  as 


Pa  ■ j'* dz  A(z)  <  E* 


E>  -  A  J°dz  A(z)o±(z)|E±(z)|2 


(19) 


The  F-region  absorption  coefficient  is  then  defined  as 

D  *  -  1  -  |Rl2  ,  (20) 

^in 

where  p^n  *  P(o)A(o)  is  the  net  HF  power  input  at  the  bottom  of  the  F-region. 
Using  equation  (18)  for  R  results  in 

rzo 

n  -  1  -  exp[-4  /dz  a  (z)]  . 

X  ± 


(21) 


According  to  equation  (11),  o±(z)  Is  a  function  of  both  the  collision 
frequencies  and  the  HF  frequency.  Therefore,  the  absorption  efficiency  n  is 
aoderately  dependent  on  temperature  and  density  variations  via  the  variation 
In  collision  frequencies,  but  strongly  dependent  on  the  variation  in  HF 
frequency,  as  shown  In  section  III  C. 

D.  Method  of  Solution 

The  transport  equations  (1)  -  (3)  together  with  the  scheme  for  evaluating 
the  heating  source  Q  explained  in  part  C  of  this  section  constitute  the  mathe¬ 
matical  formulation  of  the  HF  heating  model.  The  mathematical  structure  con¬ 
sists  of  a  set  of  coupled  nonlinear  partial  differential  equations  for  the 
electron  density  Ne(z,t)  and  electron  temperature  Te(z,t)  whose  solutions  must 
be  obtained  numerically.  It  should  be  noted  that  the  density  equation  is 
highly  sensitive  to  temperature  variations,  while  the  equation  for  the  temper¬ 
ature  is  not  strongly  affected  by  the  small  variations  in  density. 

To  formulate  a  well  posed  numerical  problem  a  set  of  boundary  conditions 
must  be  prescribed  at  z  ■  0  and  at  some  large  value  of  z  -  zu,  well  beyond 
that  corresponding  to  the  critical  density  of  the  F-region.  At  the  bottom  of 
the  F-region  (z  ■  0)  we  assume  that 

Te(z-0,t)  -  Teo(z-0)  -  Tn  .  (22) 

This  choice  is  quite  reasonable  physically,  because  at  the  lower  heights  the 
electrons  are  strongly  coupled  to  the  neutrals  and  thus  loose  their  energy 
quite  rapidly.  In  addition,  this  choice  is  consistent  with  models  for  D  and 
E-regions  heating  by  HF  waves  [Tomko,  1981].  It  is  also  assumed  that  at  the 
lower  boundary 
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Ne(z-0,t)  -  Neo(z«0)  .  (23) 

This  condition  is  justified  because  the  density  changes  produced  by  HF  heating 
are  maximized  near  the  reflection  layer  (z  ■  zQ)  and  even  there  they  tend  to 
be  small  (of  the  order  of  few  percent,  as  is  shown  in  section  III  A). 

At  the  upper  boundary  z  =  zu  we  demand  that  the  gradients  of  the 
zero-order  profiles  remain  unchanged,  i.e. , 

[='•' . • 

(24) 

[sV-4  ■  [Uw],  • 

This  choice  is  motivated  by  the  physical  idea  that  the  particle  and  heat 
fluxes  at  altitudes  between  800  km  to  1000  km  are  not  likely  to  be  influenced 
by  HF  pulses  that  reflect  around  250-300  km  altitudes.  However,  from  a  rigor¬ 
ous  point  of  view  the  upper  boundary  conditions  are  unknown  since  the  plasma 
there  tends  to  be  collisionless,  hence  it  does  not  admit  a  fluid  transport  de¬ 
scription.  In  practice,  the  upper  boundary  introduces  numerical  reflections 
which  can  invalidate  some  results,  in  particular,  the  density  changes  produced 
by  long  HF  pulses.  However,  all  the  results  reported  in  the  present  work  are 
obtained  in  the  regime  where  the  sensitivity  to  the  upper  boundary  is  negligi¬ 
ble. 

Equations  (l)-(3)  are  solved  using  a  predictor-corrector  form  of  the 
time  averaged  Crank-Nicholson  tridiagonal  method  [Douglas  and  Jones,  1963], 

The  spatial  grid  scale  is  chosen  so  that  it  resolves  heights  as  small  as  250  m, 
while  the  smallest  time  scales  resolvable  are  of  the  order  of  100  ms.  Non- 


H 


1 


> 
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linear  coefficients  are  updated  to  half  the  time  step  before  the  variables  are 
time  averaged  over  a  full  time  step.  The  new  values  of  density  and  tempera¬ 
ture  are  then  fed  back  into  all  the  dependent  parameters  of  the  analysis. 
Specifically,  the  feedback  of  the  heating  on  the  source  term  Q  shews  through 
the  collision  dependent  parameters,  i.e.,  the  conductivity  o+  and  the  absorp¬ 
tion  o+.  The  effect  of  the  density  variation  on  the  spatial  shape  of  the 
HF  field  near  the  reflection  layer  is  neglected  because  the  density  changes 
are  small.  Also,  in  our  computations  we  average  the  spatial  sinusoidal  varia¬ 
tions  resulting  from  the  interaction  between  the  incident  and  reflected  waves, 
the  reason  being  that  their  scale  length  is  usually  smaller  than  the  collision 
length  of  the  particles,  making  the  transport  analysis  insensitive  to  their 
effect.  However,  this  approximation  does  not  affect  the  magnitude  and  shape 
of  the  large  main  Airy  pattern  lobe,  especially  on  the  evanescent  high  density 
side,  which  is  important  in  calculating  the  density  modification. 
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III.  NUMERICAL  RESULTS 

A.  Response  to  HF  Heating 

The  HF  transmitter  characteristics  chosen  In  the  numerical  study  of  equa¬ 
tions  (l)-(3)  correspond  roughly  to  those  available  at  the  experimental  facil¬ 
ities  in  Norway  and  Alaska.  The  HF  frequency  is  taken  as  ui/2ir  >5.1  MHz,  a 
radiated  power  of  W  *  1  MW  is  chosen,  and  an  antenna  gain  G  -  20  dB  is  as¬ 
sumed.  A  pure  left  hand  circularly  polar.'  •id  wave  (corresponding  to  x-mode 
for  oblique  case)  is  considered  in  this  geometry  (VNeo  x  Bq  *  0).  The  results 
obtained  with  this  heating  mode  are  easily  extrapolated  to  the  case  of  a  fi¬ 
nite,  but  small,  angle  between  VNeo  and  Bq,  as  the  characteristics  of  the 
left  hand  wave  (x-mode)  change  slowly  with  small  angular  variations.  The 
right  hand  circularly  polarized  mode  (corresponding  to  o-mode  for  oblique 
case)  which  reflects  at  a  higher  altitude,  has  similar  Ohmic  absorption  char¬ 
acteristics  to  the  left  hand  mode;  however,  the  results  obtained  for  a  purely 
longitudinal  geometry  cannot  be  extrapolated  to  the  case  of  a  finite  angle  be¬ 
tween  VNeo  and  Bo*  For  such  a  polarization  one  must  properly  include  the  ef¬ 
fect  of  electrostatic  modes  appearing  near  the  layer  where  u>p(z)  *  ui  [Meltz  et 
al. ,  1974],  as  well  as  the  effect  of  Z-trace  propagation,  both  of  which  create 
new  absorption  channels  in  addition  to  Ohmic  heating.  The  D  and  E-regions 
transmission  factor  3  is  arbitrarily  taken  as  S  *  0.75,  which  corresponds  to 
weak  D  and  E-reglon  conditions,  a  situation  consistent  with  the  quiet  F-region 
profiles  shown  in  Figure  la. 

We  proceed  next  to  describe  the  results  obtained  from  heating  with  typi¬ 
cal  HF  pulses  90  seconds  long.  The  relaxation  of  the  F-region  is  observed  for 
60  seconds  after  the  HF  pulse  is  turned  off. 

The  altitude  dependence  of  the  local  Ohmic  heating  absorption  rate 
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Q  ■  (o_|E_j2)/2,  obtained  from  equation  (5),  Is  shown  In  Figure  4a.  In  the 
neighborhood  of  z  m  150  km,  the  absorption  Is  reduced  due  to  two  effects: 

1)  beam  divergence  and  consequent  ll/(D+z)2]  decrease  In  |E_|2;  2)  wave  atten¬ 
uation  due  to  electron  collisions  (which  decrease  with  altitude).  A  sharp  In¬ 
crease  In  Q  exists  at  the  reflection  layer  as  e_(z)  0  because  of  the  | e_ 

(z)|"l/2  dependence  of  |E_|2;  the  maximum  amplitude,  however,  is  correctly 
limited  by  the  Airy  function  as  given  in  equation  (17). 

Figure  4b  shows  the  altitude  dependence  of  the  fractional  change  in  elec¬ 
tron  temperature,  normalized  to  the  local  pre-heatiag temperature  Teo(z),  for 
different  times  in  the  evolution.  The  HF  pulse  is  turned  on  at  t  *  0.  For 
early  times,  t  <  1  s  (not  shown),  the  spatial  variation  of  the  temperature  is 
essentially  the  same  as  that  of  the  HF  power,  i.e.,  the  process  is  local  and 
is  thus  maximized  around  the  reflection  layer.  Electron  thermal  conduction 
rapidly  (within  1  s)  starts  diffusing  heat  from  regions  of  high  concentration 
near  the  reflection  layer  and  at  the  bottom  boundary  to  cooler  regions  along 
the  magnetic  field  lines  (t  ■  3  s).  The  temperature  profile  broadens  conse¬ 
quently  from  the  bottom  of  the  F-region  (z  »  0,  at  150  km)  to  the  reflection 
layer  (z  ■  zQ,  at  206  km),  with  a  small  foot  expanding  beyond  the  reflection 
layer  to  the  region  of  higher  thermal  conductivity  (t  ■  15  s).  Equilibrium  is 
attained  between  z  ■  0  and  z0  essentially  within  10  seconds  (equilibrium  time 
being  defined  as  the  time  for  the  variation  to  reach  over  95%  of  its  value), 
and  a  maximum  local  change  of  about  50%  is  produced.  After  90  s,  the  tempera¬ 
ture  of  the  whole  F-region  is  well  into  the  steady  state  equilibrium  regime. 
After  the  HF  pulse  is  turned  off,  thermal  conduction  relaxes  the  temperature 
profile  between  z  ■  0  and  zQ  essentially  back  to  its  unperturbed  level  within 
10  seconds  (t  -  105  s).  One  minute  after  the  heating  is  turned  off,  at 
t  ■  150  s,  there  is  almost  no  trace  of  HF  heating  in  the  electron  tempera¬ 


ture. 
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The  relatively  fast  relaxation  of  the  temperature  is  to  be  contrasted 
with  the  behavior  of  the  density  modification,  which  responds  on  a  much  longer 
time  scale  because  it  is  regulated  by  ambipolar  diffusion.  The  density  chan¬ 
ges  are  triggered  by  local  increases  in  plasma  pressure  resulting  from  the 
temperature  increases.  Specifically,  the  rate  of  density  change  is  governed 
by  the  second  spatial  derivative  of  the  electron  temperature.  Consequently, 
for  early  times  after  turning  on  the  HF  pulse,  the  density  profile  experiences 
an  impulse-like  driving  force  having  a  spatial  dependence  sharply  peaked  at 
the  reflection  layer  (on  the  higher  density  side).  This  effect  produces  a 
dent  in  the  density  profile  as  seen  in  Figure  4c.  This  figure  illustrates  the 
time  evolution  of  the  fractional  density  change,  normalized  to  the  local  un¬ 
perturbed  value  Neo(z),  and  it  shows  the  slow  development  of  a  "thermal  cavi- 
ton"  having  a  spatial  scale  initially  determined  by  the  width  of  the  electro¬ 
magnetic  Airy  pattern.  This  density  cavity  is  to  be  distinguished  from  other 
analogous  effects  previously  discussed  in  the  plasma  physics  literature,  e.g. , 
ponderomotive  density  cavities  which  can  be  generated  at  short  scales  by  elec¬ 
trostatic  waves,  or  at  long  scale  lengths  by  standing  electromagnetic  waves. 
Actually,  it  can  be  shown  that  the  ratio  of  the  ponderomotive  force  produced 
by  the  Airy  pattern  of  (17)  to  the  kinetic  pressure  force  producing  these 
thermal  cavitons  is  of  the  order  of  (nig/mj).  As  seen  in  Figure  4c,  the  grtwth 
of  the  density  cavity  continues  throughout  the  length  of  the  HF  pulse;  eventu¬ 
ally  the  cavity  attains  a  depth  of  about  2X  of  the  unperturbed  density  and  ex¬ 
tends  over  a  length  of  ~10  km.  One  minute  after  turning  off  the  HF,  at 
t  *  150  s,  the  cavity  still  exists  near  the  cut-off  layer  at  a  level 
| 6N/Me0|  ~  IX.  The  diffusion  velocity  is  so  small  [equation  (2)]  that  the 
cavity  is  almost  static  during  the  short  time  of  the  experiment. 
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A  closer  examination  of  Figure  4c  reveals  another  density  cavity  located 
at  a  height  slightly  above  the  lower  numerical  boundary.  Although  the  depth 
of  this  lower  cavity  is  indeed  enhanced  by  our  choice  of  a  rather  abrupt  cool¬ 
ing  of  the  electron  temperature  at  the  lower  end,  the  qualitative  existence  of 
such  a  cavity  is  expected  on  first  principles.  This  second  cavity  arises  be¬ 
cause  a  second  order  spatial  derivative  in  temperature  must  develop  in  the 
transition  from  E-region  to  the  lower  F-region,  hence  an  effect  analogous  to 
the  impulsive  drive  encountered  at  the  HF  cut-off  must  also  occur  at  a  lower 
height.  A  density  "bump”  is  created  between  the  two  cavities.  It  attains  a 
level  $N/Neo  of  the  order  of  0.5%  after  90  s  of  heating  and  relaxes  at  a  rate 
much  slower  than  the  cavities. 

Figure  4d  shows  the  electron-neutral  and  electron-ion  cooling  losses  be¬ 
fore  HF  heating  (t  ■  0),  after  HF  heating  is  on  (t  *  90  s),  and  one  minute  af¬ 
ter  HF  heating  is  turned  off  (t  ■  150  s).  It  is  interesting  to  note  the  large 
temperature  dependence  of  the  inelastic  electron-neutral  cooling  processes. 
Their  presence  plays  an  important  role  in  limiting  the  heating  of  the  elec¬ 
trons.  This  effect,  together  with  thermal  conduction,  limits  the  production 
of  the  thermal  runaway  process  discussed  by  Perkins  and  Roble  [1978], 

Figure  4e  shows  the  behavior  of  the  transverse  d.c.  Pedersen  electrical 
conductivity  op(z,Te)  and  the  parallel  d.c.  electrical  conductivity  0|(z,Te) 
for  the  electrons  before  HF  heating  (t  *  0),  after  heating  is  on  (t  *  90  s), 
and  one  minute  after  it  is  turned  off  (t  =  150  s).  We  do  not  show  the  Hall 
conductivity  because  it  is  not  sensitive  to  temperature  changes,  and  remains 
essentially  constant  in  the  F-reglon  because  the  density  changes  are  small. 
Equations  (10)  show  that  when  u>  ■  0  (d.c.  case),  cp  «  vg  and  Oj  *  ve~*»  yield- 


|  6o.(z,T  )/c.  (z)  |  «  ,V(.Z) . .  |  5a  (z,T  )/o  (z)|  , 

1  e  lo  ve(  z  ,Te)  P  e  po 
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where  0^(2),  0|o(z),  veo(z)  are  the  initial  values  before  heating.  Heating 
decreases  ve  in  the  middle  and  upper  F-reglon  where  electron-ion  collisions 
are  dominant.  It  Increases  ve  in  the  lower  F-region  where  electron-nitrogen 
momentum  transfer  collisions  are  dominant.  This  behavior  is  observable  in  the 
variation  of  Op  and  0|  in  Figure  4e.  It  is  of  interest  to  estimate  the  total 
change  in  the  electrical  resistance  along  a  magnetic  flux  tube  of  fixed  cross 
section  in  the  F-region.  Alternatively,  the  variation  in  the  height- 
integrated  parallel  resistivity  is  given  by: 

Ar  r  Zu  _  /  '  r  zu 

Ji-  r_L _ -_L_id,/  (26) 

rlo  0  L  0|(z.Te)  J  /  ■'o  °lo(z> 


where  rj  is  the  parallel  resistance  of  the  flux  tube  and  r|c  is  its  value  be¬ 
fore  heating  (the  effect  of  the  ions  being  negligible).  We  find  that  after  90 
s  of  heating,  the  parallel  resistance  in  the  F-region  has  decreased  by  about 
4Z.  The  corresponding  variation  in  the  height  integrated  Pedersen  resistivity 
is  negligible  as  it  is  dominated  by  the  unchanged  ion  component.  The  net  ef¬ 
fect  may  alter  the  flow  of  naturally  occurlng  auroral  currents  (Birkland  cur¬ 
rents)  by  favoring  their  parallel  flow  in  the  F-region  and  lowering  the  alti¬ 
tude  of  their  return  path  at  the  E-region  heights.  Such  interactions  are  of 
physical  interest  and  warrant  further  study. 

A  common  feature  in  the  results  shown  in  Figure  4  is  that  maximal  varia¬ 
tion  occurs  near  the  reflection  layer,  where  maximum  Ohmic  absorption  takes 
place.  In  Figure  5  we  show  the  time  variation  at  the  reflection  layer  zQ 
(corresponding  to  an  altitude  of  206  km)  of  the  temperature,  the  density,  and 
the  electron  Pedersen  conductivity,  respectively.  The  rapid  rise  and  decay 
time  of  the  temperature,  as  previously  mentioned,  is  more  visible  in  Figure  5a 
and  is  found  to  possess  an  e-foldlng  time  of  about  1  s.  The  slow  diffusic  of 


the  density  is  shown  in  Figure  5b.  Contrary  to  the  temperature,  no  steady 
state  or  return  to  background  conditions  Is  observed,  only  slow  rise  and  de¬ 
cay.  The  variation  of  the  electron  Pedersen  conductivity  follows  closely  that 
of  the  temperature  as  seen  in  Figure  5c.  The  long  term  steady  state  is 
altered  only  by  the  very  small,  continuous  variation  of  the  density.  The  cor¬ 
responding  variation  of  the  parallel  conductivity  can  be  inferred  from  Equa¬ 
tion  (25),  whence  it  is  expected  to  look  similar  to  that  of  the  Pedersen  com¬ 
ponent,  but  positive  in  magnitude. 


B.  Power  Dependence 

The  numerical  results  discussed  thus  far  correspond  to  the  effects  pro¬ 
duced  by  a  fixed  average  radiated  power  W  *  1  MW.  Due  to  the  nonlinear  depen¬ 
dence  of  the  various  collisional  processes  on  electron  temperature,  the  result 
ing  values  of  the  temperature  variation,  depth  of  density  cavities,  and  other 
dependent  parameters  are,  in  general,  nonlinear  functions  of  the  radiated  pow¬ 
er.  To  study  such  nonlinear  behavior  we  proceed  to  vary  the  radiated  power 
with  parameters  similar  to  those  used  in  section  III  A,  specifically  a  fixed 
frequency  of  5.1  MHz  and  a  pulse  length  of  90  s.  The  time  evolution  of  the 
temperature  and  density  at  the  reflection  layer  zQ  (corresponding  to  an  alti¬ 
tude  of  206  km)  are  shown  in  Figure  6,  and  the  dependence  of  these  variations 
on  the  radiated  power,  after  90  s  of  heating,  is  shown  in  Figure  7.  The  vari¬ 
ation  of  the  electron  Pedersen  conductivity  is  also  shown  in  Figure  7  (curve 
c).  Besides  the  larger  variations  achieved  with  higher  power.  Figure  6  also 
shows  that  the  rise  and  decay  equilibrium  times  for  the  temperature  are  power 
dependent,  the  higher  the  power  the  shorter  the  time  (~10  s  for  W  =  1  MW,  and 
~5  s  for  W  «  5  MW).  This  effect  is  stronger  as  one  moves  toward  lower  F- 
region  altitudes  (not  shown).  It  is  seen  from  Figure  7  that  clear  departure 
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from  linearity  sets  in  at  powers  larger  than  W  =  1  MW.  Saturation  of  both  the 
temperature  and  the  conductivity  (curves  a  and  c)  seem  to  be  attainable  at 
power  levels  higher  than  W  =  10  MW.  Curve  b  shows  that  density  cavities  with 
depths  of  about  4%  of  the  background  density  may  be  attained  at  power  levels 
W  >  5  MW. 

The  spatial  variation  of  the  electron  temperature  for  different  power 
levels  is  shown  in  Figure  8  where  it  is  seen  that  temperature  variations  over 
120%  are  attained  in  the  lower  F-region  for  W  >  5  MW.  Also,  heat  conduction 
for  W  >  5  MW  can  increase  the  temperature  of  the  upper  F-region  by  almost  10% 
of  its  background  value. 

C.  Frequency  Dependence 

The  next  variation  we  investigate  is  the  dependence  of  the  ionospheric 
modifications  on  the  choice  of  HF  frequency  for  fixed  average  radiated  power 
W  *  1  MW,  and  with  other  parameters  similar  to  those  used  in  section  III  A. 

All  cases  presented  correspond  to  overdense  heating  (where  (o  <  (^max)*  Fur¬ 
thermore,  to  illustrate  the  degree  of  variability  which  may  be  encountered  in 
field  experiments,  we  consider  two  model  ionospheres  having  identical  height 
dependence  (self-similar  profiles),  but  with  different  critical  frequencies 
a)P2max*  The  first  model,  a  (shown  in  Figure  la),  has  Neo(0)  =  5  x  10^cra~\ 
and  the  second  model,  b,  has  Neo(0)  =  lO^cm-^. 

Figure  9  shows  the  dependence  of  the  F-region  absorption  coefficient 
n  on  HF  frequency  for  the  two  model  ionospheres  a  and  b.  The  solid  curves 
correspond  to  conditions  prior  to  application  of  the  HF  pulse  while  the  dashed 
curves  represent,  the  situation  aftei  90  s  of  HF  heating.  The  noticeable  depen¬ 
dence  In  Figure  9  is  the  increase  of  the  absorption  with  HF  frequency,  a  re¬ 
sult  associated  with  the  longer  absorption  path  length  travelled  by  the  higher 
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frequency  pulses.  Depending  on  the  frequency  and  the  density  profile,  the  ab- 
sorption  coefficient  n  varies  between  20-60%  for  the  quiet  conditions  of  our 
model.  For  neutral  temperatures  corresponding  to  higher  solar  activity  the 
upper  bound  on  n  may  be  over  80%.  Another  dependence  seen  in  Figure  9  is  that 
after  90  s  of  heating  the  absorption  coefficient  drops  by  about  5-10%  from  its 
initial  value.  This  behavior  is  to  be  contrasted  with  that  occuring  for  fre¬ 
quencies  that  reflect  near  the  lower  F-region  where  the  absorption  increases 
after  the  application  of  the  HF,  a  characteristic  of  the  D  and  E-regions. 

The  relative  variation  of  the  temperature  produced  at  the  reflection 
layer  after  90  s  of  HF  heating  (steady  state)  is  plotted  against  frequency  in 
Figure  10  for  the  two  model  ionospheres  a  and  b.  The  behavior  is  character¬ 
ized  by  temperature  increases  between  50-90%  at  lower  frequencies,  which  cor¬ 
respond  to  reflection  in  the  highly  collisional  lower  F-region.  The  tempera¬ 
ture  variation  decreases  as  the  HF  frequency  increases  and  reaches  a  nearly 
constant  value  of  [ 6Te/Teo(z) ]  ~  30%,  for  both  ionospheres  a  ard  b.  This 
flattening  effect  arises  for  HF  frequencies  which  reflect  in  the  region  loca¬ 
ted  below  the  F-region  peak  where  the  total  electron  collision  frequency 
( ve  =  ve£  +  ven)  exhibits  a  weak  spatial  dependence. 

Figure  11  shows  the  frequency  dependence  of  the  peak  relative  change  in 
density  associated  with  the  cavities  generated  at  the  reflection  layer  after 
90  s  of  heating.  The  overall  dependence  of  the  density  cavity  follows  essen¬ 
tially  the  same  behavior  exhibited  by  the  temperature  as  expected,  since  their 
formation  is  triggered  by  the  sharp  temperature  gradients  at  the  reflection 
layer.  An  additional  feature  is  observed  in  Figure  11.  Instead  of  approach¬ 
ing  a  flat  frequency  dependence  at  the  higher  frequencies  as  in  Figure  10,  the 
magnitude  of  [ 6 Ne / Ng 0 ( z ) ]  tends  to  increase.  This  secondary  effect  is  associ- 
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ated  with  frequencies  that  reflect  at  higher  altitudes  where  the  exponentially 
decreasing  ion-neutral  collision  frequency  is  responsible  for  increasing  the 
rate  of  ambipolar  expansion. 

Finally,  a  general  observation  regarding  the  previous  results.  The 
curves  in  Figures  9,  10,  and  11  for  ionosphere  b  are  almost  similar  to  those 
obtained  for  ionosphere  a  but  shifted  towards  higher  frequencies.  This  is  an 
indication  that  the  depth  of  penetration  inside  the  F-region,  rather  than  the 
frequency  itself,  is  the  major  factor  in  determining  the  magnitude  of  the  var¬ 
iation  of  the  parameters. 

D.  HF  Modulation  Pulses 

The  effect  of  modulating  the  F-region  with  HF  pulses  repeated  at  low  fre¬ 
quency  (ULF  range)  is  a  topic  of  general  interest  which  we  investigate  here 
using  a  1  Hz  square  amplitude  modulation  of  the  HF  pump.  All  HF  parameters 
are  similar  to  those  used  in  section  III  A,  specifically,  the  frequency  is  set 
at  5.1  MHz  and  the  average  radiated  power  is  initially  assumed  to  be  W  =  1  MW. 
Figure  12a  shows  the  modulation  of  the  electron  temperature  and  Figure  12b  the 
corresponding  modulation  of  the  electron  Pedersen  conductivity,  from  which  the 
modulation  of  the  parallel  conductivity  can  also  be  inferred  using  equation 
(25).  With  the  power  held  at  W  *  1  MW  we  modulate  from  steady  state  condi¬ 
tions  after  30  s  of  HF  heating.  We  obtain  about  2.5%  modulation  of  the  con¬ 
ductivity  corresponding  to  5%  modulation  of  the  temperature.  If  we  reduce  the 
pre-modulation  heating  time  to  5  s  keeping  W  =  1  MW  and  then  modulate  at  1  Hz 
(i.e.  ,  we  do  not  modulate  from  steady  state),  the  result  is  the  same.  If  we 
double  the  power  to  W  *  2  MW,  again  with  5  seconds  heating  before  the  1  Hz 
modulation,  we  notice  that  the  initial  variation  of  both  the  temperature  and 
conductivity  is  larger,  however  the  conductivity  modulation  attained  is  only 


augmented  to  approximately  3%  corresponding  to  about  9%  temperature  variation 
The  modulations  are  evaluated  at  the  reflection  layer  of  the  HF  wave,  hence 
they  correspond  to  maximal  variation  as  is  shown  in  section  III  A.  It  is  sug 
gestive  then  that  in  order  to  enhance  the  modulation  of  the  conductivity  one 
should  use  a  pump  frequency  that  reflects  in  the  lower  F-region  where  higher 
electron  collision  rates  prevail.  This  idea,  however,  is  only  partially  cor¬ 
rect.  As  may  be  seen  from  Figure  4e,  there  is  a  location  in  the  F-region 
where  the  Pedersen  conductivity  increases  (oy  decreases)  with  temperature  in¬ 
stead  of  decreasing  (ay  increasing)  as  is  the  case  in  the  rest  of  the  F- 
region.  This  behavior  is  related  to  changing  from  a  domain  where  electron- 
nitrogen  collisions  dominate  to  one  where  electron-ion  collisions  dominate. 
Therefore,  there  exists  a  layer  in  the  lower  F-region,  between  160-200  km 
altitude  for  the  present  case,  where  both  effects  tend  to  cancel  each  other 
and  in  the  neighborhood  of  which  large  temperature  variations  produce  no  sig¬ 
nificant  modulation  of  the  conductivity. 

The  effect  of  the  modulation  on  the  density  Is  to  keep  the  depth  of  the 
density  cavities  nearly  constant,  i.e.  ,  at  the  level  attained  just  before  mod¬ 


ulation  starts 


IV  CONCLUSIONS 


We  have  presented  a  model  for  the  Ohmic  heating  by  HF  waves  of  the  polar 
Ionosphere  F-reglon.  The  geometry  of  this  high-latitude  allows  a  self- 
consistent  analysis  of  overdense  heating.  The  ability  to  analytically  model 
the  spatial  behavior  of  the  HF  electric  field  proves  essential  in  obtaining 
accurate  estimates  of  the  changes  produced  in  the  temperature  and  density  pro¬ 
files. 

Ohmic  absorption  varying  from  202  to  802  is  attainable,  depending  mainly 
on  the  HF  frequency  and  solar  activity  cycle.  Maximum  ionospheric  modifica¬ 
tion  is  found  to  occur  near  the  reflection  layer  where  Ohmic  absorption  is 
largest.  This  behavior  is  in  agreement  with  experimental  data  [e.g. ,  Gordon 
and  Carlson,  1974],  and  results  in  the  formation  of  localized  thermal  cavi- 
tons.  Maximum  electron  temperature  increases  in  the  range  of  20-1202  are  ob¬ 
tained  for  HF  power  levels  ranging  between  0.2-5  MW,  and  exhibit  a  fast  rise 
and  decay  equilibrium  time  of  7-10  s  (depending  on  altitude)  for  powers  near 
W  ~  1  MW,  and  of  2-5  s  (depending  on  altitude)  for  higher  powers  W  >  5  MW. 
These  times  are  shorter  than  those  for  mid-latitude  heating,  however  the 
general  behavior  is  consistent  with  experimental  observations  [e.g.,  Gordon 
and  Carlson,  1974;  Mantas,  et  al. ,  1981 ]  as  well  as  with  previous  calculated 
values  [Meltz  and  Lelevier,  1 9 70 ] . 

The  thermal  density  cavities  produced  at  the  reflection  layer  are  a  few 
kilometers  in  extent  along  the  earth’s  magnetic  field  and,  depending  on  HF 
power  and  pulse  length,  can  attain  a  depth  of  1-10%  of  the  local  background 
density.  These  values  agree  with  the  density  changes  inferred  from  experi¬ 
ments  [e.g.,  Minkoff,  1974;  Duncan  and  Behnke,  1978].  The  response  time  of 
these  termal  cavitons  is  slow,  of  the  order  of  minutes.  Hence  they  remain 
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for  a  considerable  time  after  the  HF  pulse  is  turned  off.  This  behavior  re¬ 
sults  in  a  memory-type  effect,  whereby  in  a  repeated  HF  pulse  scenario  with 
off  time  of  seconds  to  few  minutes,  the  ionosphere  "remembers"  the  effect  of 
previous  pulses.  Any  pre-existing  density  cavity,  whether  produced  by  a  pre¬ 
vious  HF  pulse  or  naturally  occurring,  can  be  stimulated  into  growing  by  sub¬ 
sequent  HF  pulses.  For  instance,  the  depth  of  the  density  cavity  formed  with 
1  Hz  modulation  presented  in  section  III  0  can  be  maintained  at  essentially  a 
constant  level.  The  existence  of  such  cavities  near  the  reflection  layer  can 
in  turn  modify  the  reflection  properties  of  test  signals,  such  as  ionosondes. 

Such  an  effect  is  consistent  with  the  experimental  observation,  first  seen  in 
Platteville,  of  artificial  spread-F  persisting  long  after  the  HF  pump  has  been 
turned-off  (Utlaut  and  Violette,  1974], 

The  present  study  strongly  suggests  that  the  large  estimates  of  density 
modification  calculated  from  steady  state  analysis  of  the  temperature,  via 
the  pressure  equilibrium  condition,  are  never  reached  in  practice  [e.g. ,  Far¬ 
ley,  1963;  Gurevich,  1967;  Ignat’ev  et  al. ,  1977],  The  reason  is  that  it 
takes  several  minutes  for  pressure  equilibrium  to  be  attained.  However,  on 
these  time  scales  plasma  drifts  would  shift  the  heated  volume  out  of  the  HF 
beam  cross-section  [Holway  et  al. ,  1978],  thus  limiting  the  maximum  value  of 
the  density  changes  produced.  The  present  study  tends  to  show  a  general 
agreement  with  the  estimates  of  the  time  dependent  analysis  of  Meltz  and  Lele- 
vier  [1970]  and  Meltz,  Holway,  and  Tomljanovich  [1974]  which  are  obtained  for 
a  mid-latitude  ionospheric  model. 

It  has  also  been  found  in  the  present  study  that  HF  pulses  can  change  the 
d.c.  electrical  resistance  along  a  magnetic  flux  tube,  decreasing  it  by  about  4%. 
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This  suggests  the  possibility  that  parallel  auroral  Birkland  currents  may  be 
modified  with  HF  heating  by  forcing  their  return  path  towards  lower  altitudes. 

Finally,  a  useful  result  can  be  inferred  from  the  study  reported  here 
for  ionospheric  F-region  experiments  that  aim  to  investigate  F-region  wave- 
plasma  phenomena  such  as  wave  beating  [Shoucri  et  al. ,  1982] ,  parametric  decay 
instabilities  [Perkins  et  al. ,  1974] ,  electrostatic  plasma  line  effects  [Gor¬ 
don  and  Carlson,  1974] ,  etc.  Such  experiments  use  characteristic  pulse 
lengths  of  the  order  of  a  few  msec  to  a  few  hundred  msec.  According  to  the 
present  results  such  time  scales  are  extremely  short  compared  to  the  longer 
time  scales  for  Ohmic  heating,  cooling,  and  density  diffusing.  Consequently, 
the  HF  pump  is  expected  to  have  a  negligible  effect  on  the  background  iono¬ 
sphere  during  such  short  time  experiments.  However,  the  use  of  repetitive 
pulses  at  small  time  intervals  to  conduct  such  experiments  as  well  as  the  use 
of  pulse  lengths  of  a  few  seconds  can  produce  significant  changes  in  the  back¬ 
ground  ionosphere. 
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FIGURE  CAPTIONS 


1_  The  ionospheric  F-region  model:  a)  zero  order  Chapman  electron  den¬ 
sity  and  temperature  profiles;  b)  density  profiles  of  various  neu¬ 
tral  components;  c)  elastic  momentum  transfer  collision  profiles. 

2_  Different  electron  cooling  mechanisms  in  the  F-region  evaluated  at 
equlllbriim  prior  to  HF  heating  and  due  to:  elastic  collisions  be¬ 
tween  electron-neutrals  (e-n)R^  and  electron-ions  e-i;  inelastic 
collisions  between  electrons  and  rotational  levels  e-N2R  and  e-02R, 

vibrational  levels  e-N2^  and  e-02V»  oxygen  fine  structure  of  ground 

3P  1 

level  e-0  ,  and  oxygen  first  excited  state  e-0  D;  e-n  is  the  total 

electron-neutral  losses. 

3  Schematic  diagram  for  the  HF  heating  of  the  polar  F-region.  The 

distribution  of  the  squared  HF  electric  field  is  shown,  and  its  fine 
structure  near  the  reflection  layer  is  magnified  in  the  boxed  area. 

4^  Ionospheric  response  to  a  5.1  MHz,  left  handed  circularly  polarized 
(L-  or  x-mode)  HF  pulse  of  length  t  ■  90  s;  radiated  power  W  -  1  MW, 
antenna  gain  G  -  20  dB,  D  and  E-regions  transmission  factor 
8  ■  0.75,  after-pulse  cooling  time  60  s.  Shown  is  the  altitude  de¬ 
pendence  at  different  times  of:  a)  Ohmic  heating  Q  *  o_|e_|2/2; 
b)  electron  temperature  variation,  STe  *  Te(z,t)  -  Teo(z);  c)  plasma 
density  variation,  6Ne  -  Ne(z,t)  -  N0O(z);  d)  electron  cooling 
losses  due  to  neutrals  Len  and  ions  Le^;  e)  variation  of  the  elec¬ 
tron  Pedersen  and  parallel  conductivities,  op  and  0|  respectively. 

5_  Time  variation  at  the  reflection  layer  (zQ  =  206  km,  measured  from 
the  ground)  of:  a)  electron  temperature;  b)  plasma  density;  c)  elec¬ 
tron  Pedersen  Conductivity.  Same  HF  pump  parameters  as  in  Figure  4. 


6^  Time  variation  at  the  reflection  layer  (zQ  =  206  km,  measured  from 
the  ground)  of:  a)  electron  temperature,  and  b)  plasma  density,  for 
different  HF  radiated  powers  W  -  0.2,  1,  and  5  MW.  Other  HF  pump 
parameters  are  similar  to  Figure  4. 

7_  Power  dependence  at  the  reflection  layer  (zQ  =  206  km,  measured  from 
the  ground)  of  the  electron  temperature,  plasma  density,  and  elec¬ 
tron  Pedersen  conductivity  for  the  HF  pump  parameters  of  Figure  4. 

8_  Spatial  changes  of  the  electron  temperature  for  different  HF  power 
levels. 

9_  Frequency  dependence  of  the  Ohmic  absorption  coefficient  n  before 
(t  ■  0)  and  after  heating  (t  “  90  s)  for  two  self-similar  density 
profiles  having:  (a)  Neo(0)  -  5  x  10^  cm-3,  and  (b)  Neo(0)  -  105cm-3 
W  -  1  MW,  G  -  20  dB,  L-  or  X-mode. 

10  Dependence  of  the  electron  temperature  changes  at  the  reflection 
height  on  frequency,  after  90  s  of  heating.  Results  are  shown  for 
two  self-similar  density  profiles  having  (a)  Neo(o)  -  5  x  lO^cm-3, 
and  (b)  Neo(o)  -  105cm-3;  W  -  1  MW,  G  -  20  dB,  L-mode,  B  -  0.75. 
Temperature  values  are  normalized  by  the  pre-heating  temperature  at 
the  respective  reflection  heights. 

1 1  Dependence  of  the  plasma  density  changes  at  the  reflection  height  on 
frequency,  after  90  s  of  heating.  Results  are  shown  for  two  self¬ 
similar  density  profiles  having  (a)  Neo(0)  *  5  x  lO^cm-3,  and 

(b)  Neo(0)  -  105cm-3;  W  -  1  MW,  G  -  20  dB,  L-mode,  B  -  0.75.  Den¬ 
sity  values  are  normalized  by  the  pre-heating  densities  at  the 
respective  reflection  heights. 


Figure  12  Electron  temperature  (a)  and  conductivity  (b)  modulation  by  a  1  Hz 
pulsed  HF  wave  after  30  s  of  heating  at  W  *  1  MW,  after  5  s  of  heat¬ 
ing  at  W  -  1  MW,  and  after  5  s  of  heating  at  W  *  2  MW;  f  ■  5.1  MHz, 
G  ■  20  dB,  L-mode.  Values  calculated  at  the  reflection  layer 
(zQ  =  206  km,  measured  from  the  ground).  Note  the  change  in  time 
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